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Получены аналитические зависимости и приведены результаты расчета численных значений распределения температур по по- 
верхности низкотемпературного композиционного электрообогревателя, которые могут быть использованы при проектирова- 
нии аналогичных изделий в различных областях сельского хозяйства и промышленности. 


Широкое использование низкотемпературных 
(300. ..400 К) композиционных электрообогревате- 
лей (НТКЭ) и обогревательных систем на их осно- 
ве позволяет осуществить массовое внедрение 
энергосберегающих технологий на предприятиях 
агропромышленного комплекса: электрообогрев 
индукционных счетчиков в щитах раздельного уче- 
та электроэнергии, поверхностно-распределитель- 
ный обогрев молодняка животных и птицы, подог- 
рев зерна на зерноперерабатывающих комплексах, 
обогрев задвижек, устьев скважин и трубопроводов 
различного назначения. 

Однако, в каждом конкретном случае возникает 
проблема рационального выбора типа и конструк- 
ции НТКЭ, отвечающего широкому спектру требо- 
ваний по обеспечению теплоэлектрофизических 
характеристик, электробезопасности, экономич- 
ности, стойкости к агрессивным средам. В связи с 
этим предлагается метод расчета распределения 
температуры по поверхности НТКЭ в зависимости 
от его конструктивных параметров. 

Распределение температуры по поверхности 
электрообогревателя (рис. 1, а) является его важ- 
нейшей технической характеристикой. 


Расчёт тепловой проводимости рассматривае- 
мой системы [1] позволяет получить ряд необходи- 
мых параметров для оценки распределения темпе- 
ратуры по поверхности электрообогревателя, но не 
решает до конца поставленную задачу. 

Определение распределения температуры по 
всей поверхности электрообогревателя приводит к 
довольно большому численному массиву. Для 
практических целей обычно достаточно опреде- 
лить разность температур АТ, А Т ъ А Т ъ АТ Ь АТ 4 
между характерными точками соответственно А 2 и 
А 3 , А х и А 4 , А 5 и А.,, А 6 и А 5 , А- и А 6 (рис. 1 , б). 

В связи с этим цель настоящей работы - строгое 
решение поставленной задачи и вычисление на его 
основе значений разностей указанных температур 
в широком диапазоне изменения геометрических 
параметров системы. 

Для расчёта распределения температуры по пове- 
рхности указанной системы при стационарном про- 
цессе теплового излучения воспользуемся методи- 
кой, изложенной в работе [2], основанной на непос- 
редственном определении напряжённости поля и 
конформных преобразований с учётом математичес- 
кой аналогии между потенциальными полями. 
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а) б) в) 


Рис. 1. Плоскопараллельная система: а) исходная система тел; 6) расчетная модель; в) система тел в отображенной плоскос- 
ти; 1 — горячее тело; 2 — холодная поверхность тела; 3 — теплоизолирующий материал 


Примем, как и в работе [ 1 ] , температуру горяче- 
го тела Т ь холодного Т ъ разность температур меж- 
ду ними А Т, коэффициент теплопроводности изо- 
лирующего материала А, а коэффициент теплопро- 
водности вне изолирующего материала Я,=0. 

По оси х (рис. 1 , б) введём непроницаемую плос- 
кость, примем заштрихованную часть расчётной 
модели за часть плоскости комплексного перемен- 
ного 2 и конформного, отобразим на неё верхнюю 
полуплоскость нового комплексного переменного 
С, (рис. 1, в) при соответствии точек исходной и 
отображённой плоскостей. Система уравнений, 
связывающая конструктивные размеры расчётной 
модели с координатами точек а-а 4 в отображённой 
плоскости, полученной в работе [1] с помощью ин- 
теграла Кристоффеля-Шварца, имеет вид 


В уравнениях (4-7) интегралы являются гипе- 
рэллиптическими и должны находиться численны- 
ми методами. Уравнениями (1-3) решена геомет- 
рическая часть задачи. Рассмотрим физическую 
сторону задачи. Модуль потока тепла в точках ^ 
плоскости расположения пластин к гП (^) и раз- 
ность температур АТ между первым и вторым тела- 
ми определяют по формулам [1-3]: 

Ѵ.«>- Л - 


%-а; К-аЖ-а 4 


“к + 1 “ 3 

| \__ й (Ж=А\- 


I % - а 1 ; - «з б; - а 4 


о. / 1 + / 3 
/ / 1 + / 3 + / 5 

ь _ Л . 

/ 7 ] + І 3 + / 5 


где а ь а' к - координаты краев ^-ого тела (пласти- 
ны); А - постоянная; п - число электродов. 

Из ур. (9) следует 

АТ 


7 і +/ з = !і 


I («Г -Г) 
(аІ-^)(аІ-е)(аІ-^ 2 ) 


(4) 


/ = [ | 

4 |Ѵ(а 2 2 -^)(а(-Г)(а 4 2 -Г) 

/ А I — К 

5 


(Г-^ 2 ) 

І(Г-а 2 )(^ 2 -а 2 )(а 2 -Г) С 


І ^ 2 -а 2 2 )(а 3 -^)(а 4 -^) 

Так как при конформных преобразованиях тем- 
пература в сходственных точках остаётся неизмен- 
ной, то разность температур А Гмежду точками Л' 2 и 
А\ (рис. 1, в) соответствует разности температур 
между точками А 2 и А 3 (рис. 1,6). 

Аналогично с учётом выражения (10) относи- 
тельные разности температур находят по формулам 
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где 


А = 


ЛА = А. 

АГ 1 ’ 
АТ I ’ 

ЛА = А> 

АГ / ’ 

4 


7 в = 1 

а 4 

а 4 

7 » = / 


А„ = 


= 1 


^(% 2 -а 2 _)(%-а 3 )(%-а 4 ) 

УІ(% 2 -а 2 )(%+а 3 )(Е; + а 4 ) 

. 

І ^ 2 ~а 2 )(^ + а 3 )(^ + а 4 ) 

СІЕ, 


, 2Л І(% 2 -а 2 2 )(%+а 3 )(%+а 4 ) 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


Интегралы, определяемые выражениями (11), 
(16) и (19), могут быть приведены к нормальной 
форме и записаны в виде 

/= 2 


А = 


А о = 


^/(«4 -а 2 )(а 2 +а 3 ) 
2 

^(а 4 —а 2 )(а 2 +а 3 ) 


—К (к) ; 

) 

(20) 

= Р{ср,к) ; 

і 

(21) 

■70ДД) , 

(22) 


где 7Г(&), Р((р,к) и Р((р,,7) - полный и неполные эл- 
липтические интегралы первого рода с модулем 


к = 


(а 4 + а 2 )(а 3 -а 2 ) 
(а 4 -а 1 )(а 3 + а 2 ) 


и аргументами 


(р = агсзш 


СІ2 “I" СІ д 


(23) 


(24) 


(р 3 = агсзіп 


(а 4 - а 2 ) ( а 3 - а 2 ) 
(а 4 + сі 2 )(а 3 +а 2 ) 


(25) 


С учётом соотношений (20-22) уравнения 
(12-15) примут вид: 


А7 1 _ Р(ср,к) 
АТ К (к) ’ 


(26) 

А Т 2 _ Ѵк - °2 ) ( а 2 + а з ) 
А7 ~~ 2 

А . 
ад’ 

(27) 

АТ 3 _ ^/(«4 ~а 2 )(а 2 +а 3 ) 
А Т~ 2 

А . 

К (к)’ 

(28) 

А7 1 _ 7(<р,Д) 


(29) 

АТ К (к) 



Определив на основании ур. (1-3) по заданным 
конструктивным размерам координаты точек а х -а 4 
в отображённой плоскости, а также содержащие в 
ур. (26—29) интегралы и выполнив необходимые 
вычисления в указанных выражениях, найдем от- 
носительные изменения температур А 7/Д 7... 
...АТ, /АТ между характерными точками электро- 
обогревателя. 

Рассмотрим некоторые параметры системы при 
предельных значениях размера к-Ь и к=Ь (рис. 2). 
Если к-Ь=о о, то 04=00 и, как следует из выражений 
(23) и (24), к= 1, ср= 0, следовательно Р(ср,к)=0, 
Щ ) =оо и согласно ур. (26) А 7/Д 7=0. Аналогично 
при этих условиях и А 77/А 7=0. 

Из приведённых данных следует, что при доста- 
точно большом расстоянии к-Ь температура точки 
А, стремится к температуре Т 2 поверхности А,А 4 хо- 
лодного тела. 

Если к-Ь= 0, то расчетную модель системы 
(рис. 2, а ) можно рассматривать в виде двух незави- 
симых частей - правой и левой. 

Как и ранее, примем заштрихованные части 
расчётной модели за части плоскостей комплекс- 


ах ІУ 



Рис. 2. Система при Ь=Ь: а) расчётная модель ; 6) система тел в отображённой плоскости для левой части системы ; в) система 
тел в отображённой плоскости для правой части системы: 1 — горячее тело ; 2 — холодная поверхность тела ; 3 — тепло- 
изолирующий материал 
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ного переменного 2 (рис. 2, а) и конформно отоб- 
разим на них верхние полуплоскости новых комп- 
лексных переменных С, (рис. 2 , б, в) при соответ- 
ствии точек исходных и отображённых плоскостей. 

Для расчётной модели (рис. 2, а) обозначения 
разности температур АТ, А Г, и А Т. сохраняются, 
разности температур между точками А, и А 4 , А % и А 5 
обозначим через А 7)’ и А Т 2 '. 

Проанализируем сначала разности температур 
для правой части системы. 

В соответствии с общими формулами ( 8 ) и (9) 
модуль потока тепла в точках Е плоскости С, 
(рис, 2 , в) и разности температур А Ги А Г/ находят- 
ся следующим образом: 

, . = : (30) 


дг = Д 


= А 


I 


5 -а х \\Е -а 2 1 
сіЕ, 

сіЕ, 


и \ 

I 




,7(ог -«’)(»: -«■) 

сіЕ, 

р;-^ 2 )(а 2 2 -Е 2 ) 

сіЕ, 


АТ, = А 


\І(а 2 ~% 2 )(а 2 ~% 2 ) ' 

сіЕ, 


(31) 


■1 




г 

сіЕ, 

4(Е 2 -аЖ-сТ) 

сіЕ, 


<№ 2 -аЖ-а \ ) 


(32) 


Приведём в выражениях (31) и (32) интегралы к 
нормальной форме 

2 /) 

АТ = АТ[ = — 7ЭД, (33) 

а , 

где ДА:) - полный эллиптический интеграл перво- 
го рода с модулем к=а 1 /а 2 . 

Поскольку слагаемые в формулах (31) и (32) 
равны между собой, с учётом последнего равенства 
легко определить, что разности температур на 
участках А г А п , Д 0 -Д 3 , А Г А Х , А.,-А 4 равны А Т/2. Это 
следует и из физических представлений при рас- 
смотрении указанной системы. 

Аналогично для левой части системы: 

\Ж= А 


л/Н - « 1 1 - «2 | 


(34) 


АТ,=А 


“2 

I 


сіе 


I 


сіЕ 


л ](а І -Е)(а 1 -Е) {^Е-аА(Е-с 2 ) 

і 

АГ 3 = А | 


сіЕ 


ат 4 = а 


% І(а 1 -Е)(а 2 -Е) ’ 
сіЕ 

\І(а 4 ~Е )( а 2 ~Е) 


(35) 

(36) 

(37) 


Так как разность температур А Т между точками 
А 7 и Д ' 8 (Д 7 и Д 8 ) равна нулю, то 

ДГ=ДГ 2 ’+ДГз+ДГ 4 =0. (38) 

Поскольку интегралы, входящие в ур. (35-37), 
не равны нулю, то из формулы (38) с учётом выра- 
жений (35-37) следует, что А=0 и, следовательно, 
АТ\=АТ 2 =АТ 4 =0. (39) 

Из равенства (39) следует, что при к-Ъ= 0 темпе- 
ратура на участках А Г А 6 , А^-А 5 , А 5 -А 8 равна температу- 
ре Г, резистивного материала 1 электрообогревателя. 

Очевидно, что при малом отношении (к-Ъ)/1 на 
участке А х -А ь будет также приблизительно сохранять- 
ся температура 7], а резкое уменьшение температуры 
от 7] до Т будет происходить на участке А.,-А 4 . Ука- 
занные выводы подтверждаются результатами расчё- 
тов, выполненных на ПЭВМ на языке программиро- 
вания Паскаль и представленных в табл. 1-4. 


Таблица 1. Значения АТ/АТ участков А к -А 4 и А 2 -А 3 


И/1 

аН 

Значения АТі/АТ при Ь/к, равном 



0,1 

0,2 

0.3 

0,4 

0.5 


0,1 

0,588976 

0,653322 

0.709148 

0.761292 

0,811569 


0,2 

0,642045 

0.704756 

0,752007 

0,801748 

0,849947 

1,0 

0,3 

0,680098 

0.733742 

0,783534 

0.831504 

0,877786 


0,4 

0,708818 

0.759575 

0.807669 

0.854221 

0,898654 


0,5 

0,730746 

0.779433 

0.826308 

0,871684 

0,914381 


0,1 

0,426199 

0.497386 

0.563633 

0.629730 

0,697898 


0,2 

0,467915 

0,535528 

0.600624 

0.666669 

0,734647 

1,5 

0,3 

0,499265 

0.564506 

0,628949 

0,695098 

0,763313 


0,4 

0,523667 

0,587199 

0.651295 

0,717645 

0,786002 


0,5 

0,542470 

0.604967 

0.668957 

0.735558 

0,803975 


0,1 

0,307310 

0.374826 

0,442363 

0,514376 

0,592529 


0,2 

0,337464 

0,403742 

0.472033 

0.545819 

0,626183 

2,0 

0,3 

0,360676 

0,426178 

0,495230 

0.570520 

0,652644 


0,4 

0,379022 

0,444001 

0.513820 

0.590437 

0,673991 


0,5 

0,393431 

0.458103 

0.528690 

0,606489 

0,691198 


0,1 

0,220529 

0.279872 

0.343498 

0.415784 

0,498868 


0,2 

0,241577 

0.301060 

0.366569 

0.441766 

0,528325 

2,5 

0,3 

0,257965 

0.317703 

0.384837 

0.462444 

0,551779 


0,4 

0,271010 

0,331040 

0.399640 

0,479301 

0,570901 


0,5 

0,281313 

0,341683 

0,411603 

0,492997 

0,586437 


0,1 

0,157511 

0,207473 

0,264600 

0,333458 

0,417018 


0,2 

0,171985 

0,222775 

0,282288 

0,354494 

0,442356 

3,0 

0,3 

0,183313 

0.234884 

0,296419 

0,371552 

0,462690 


0,4 

0,192391 

0.244668 

0.307959 

0.385483 

0,479365 


0,5 

0,199553 

0.252522 

0.317353 

0.396743 

0,493914 


Например, при к/1= 1, а/ 1=0, 5 и Ь/к=(к-Ь)/І= 0,5 
значения А 7,/Д 7=0, 914381; АТ^/Л 7’=0,073973; 
А Г 3 /Д 7=0, 008578; А Г 4 /А 7=0,005 106. 

В практических конструкциях электрообогре- 
вателей наиболее важно выбрать такие его 
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конструктивные размеры, при которых относи- 
тельные разности температур А Г 3 /А Ги А Г 4 / А Т бы- 
ли бы малы. Так как для электрообогревателей 
обычно выбираются отношения а/1>0,3; Ь/1> 0,5, то 
из данных табл. 3 и 4 нетрудно определить, что при 
указанных условиях (Н-Ь)/1< 0,5 значения 
ДГ 3 /ДГ<1,23 %, а ДГ 4 /ДГ<1,24 %; если же а/1> 0,3; 
Ь/1> 1 , а (Н-Ь)/1< 1 значения АГ,/АГ<15,6 %, а 
АГ 4 /АГ<1,8 %, причём указанное изменение 
А Г, /А Г происходит в основном на участке И-Ь, 
примыкающем к точке А 5 . 

Приведённые данные позволяют проектиро- 
вать электрообогреватели с заданным распределе- 
нием температуры по его рабочей поверхности. 


Таблица 2. Значения А Г 2 /Д 7 " участков и А г -А 2 


И/1 

а /1 

Значения АТ 2 /АТ при Ъ/Н , равном 



0,1 

0,2 

0,3 

0,4 

0,5 


0,1 

0,141498 

0,148332 

0,147876 

0,140662 

0,126513 


0,2 

0,150065 

0,151462 

0,144837 

0,131701 

0,111697 

1,0 

0,3 

0,154120 

0,150257 

0,139280 

0,121528 

0,097455 


0,4 

0,155464 

0,147515 

0,132679 

0,111432 

0,084733 


0,5 

0,155146 

0,143842 

0.125972 

0,102201 

0,073973 


0,1 

0,138503 

0,155622 

0,166000 

0,169538 

0,164587 


0,2 

0,150472 

0,163884 

0,170383 

0,169274 

0,158693 

1,5 

0,3 

0,158839 

0,169073 

0,172160 

0,167003 

0,151589 


0.4 

0,164818 

0,172248 

0,172347 

0,163709 

0,144276 


0,5 

0,169001 

0,174016 

0,171578 

0,160045 

0,137374 


0,1 

0,113871 

0,135425 

0.153226 

0,166907 

0,174053 


0,2 

0,124423 

0,144215 

0,160225 

0,171379 

0,174595 

2,0 

0,3 

0,132338 

0,150602 

0.164984 

0,173825 

0,173503 


0,4 

0,138413 

0,155317 

0,168241 

0,175013 

0,171558 


0,5 

0,143023 

0,158748 

0,170424 

0,175415 

0,169263 


0,1 

0,086708 

0,108248 

0,129101 

0,149171 

0,166135 


0,2 

0,094737 

0,115714 

0,136192 

0,155441 

0,170383 

2,5 

0,3 

0,100914 

0,121405 

0,141495 

0,159907 

0,172902 


0.4 

0,105766 

0,125823 

0,145550 

0,163165 

0,174343 


0,5 

0,109537 

0,129230 

0,148649 

0,165542 

0,175091 


0,1 

0,063654 

0,082977 

0,103831 

0,126665 

0,149918 


0,2 

0,069406 

0,088785 

0,110034 

0,133083 

0,155817 

3,0 

0,3 

0,073881 

0,093313 

0,1 14854 

0,138045 

0,160071 


0.4 

0,077442 

0,096915 

0,118690 

0,141908 

0,163222 


0,5 

0,080230 

0,099762 

0,121740 

0,144938 

0,165945 


Выполненные расчеты были использованы при 
проектировании и производстве НТКЭ на основе 
бутилкаучука и технического углерода промыш- 
ленных марок. Система обогрева трехфазных счет- 
чиков в щитах раздельного учета электроэнергии 
на основе разработанных электрообогревателей 
внедрена в ОАО "Первомайское Агропромэнерго" 
Алтайского края и при участии регионального 
центра ресурсосбережения Томского политехни- 
ческого университета на предприятиях агропро- 
мышленного комплекса Томской области. 
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Таблица 3. Значения АТ^/АТ участков А 6 -Д и А 2 -А 3 


Ы 

а /1 

Значения АТ 3 /АТ при Ъ/Н, равном 



0,1 

0,2 

0.3 

0,4 

0,5 і 


0,1 

0,039252 

0,037943 

0.033309 

0.026516 

0,018765 


0,2 

0,042680 

0,039802 

0.032145 

0.023680 

0,015098 

1,0 

0,3 

0,045896 

0,040026 

0.031222 

0,021402 

0,012345 


0,4 

0,049356 

0,041293 

0.030569 

0,019549 

0,010224 


0,5 

0,053322 

0,042893 

0.030187 

0,018046 

0,008578 


0,1 

0,140905 

0,144256 

0.133223 

0,111716 

0,083532 


0,2 

0,154667 

0,151035 

0.133067 

0,105667 

0,073964 

1,5 

0,3 

0,167109 

0,157189 

0.132967 

0.100536 

0,065942 


0,4 

0,179463 

0,163297 

0.133029 

0.096051 

0,059226 


0,5 

0,192470 

0,169528 

0.133174 

0,092029 

0,053544 


0,1 

0,266662 

0,283413 

0.270310 

0.234637 

0,183603 


0,2 

0,292514 

0,297732 

0,272986 

0.227055 

0,168577 

2,0 

0,3 

0,315384 

0,310148 

0,274896 

0,220202 

0,155939 


0,4 

0,337096 

0,321476 

0,276198 

0,213718 

0,144935 


0,5 

0,358376 

0,331838 

0,276760 

0.207401 

0,135237 


0,1 

0,386260 

0,420808 

0,410114 

0.365176 

0,295311 


0,2 

0,421918 

0,441224 

0.415426 

0.357258 

0,277304 

2,5 

0,3 

0,452995 

0,458341 

0.419005 

0.349602 

0,261803 


0,4 

0,481645 

0,473118 

0,421029 

0.341912 

0,247989 


0,5 

0,508584 

0,485594 

0.421454 

0,334122 

0,235750 


0,1 

0,488252 

0,540597 

0.534655 

0.485625 

0,403853 


0,2 

0,530499 

0,564698 

0.541601 

0.476495 

0,384774 

3,0 

0,3 

0,566876 

0,584359 

0.546028 

0.470164 

0,368087 


0,4 

0,599641 

0,600585 

0.548223 

0.462009 

0,353130 


0,5 

0,629553 

0,613457 

0.548215 

0.453843 

0,339127 


Таблица 4. Значения АТ^АТ участков А Г А 6 и А г -А ъ 


Н /1 

аН 

Значения Д7УДГ при Ъ/Н, равном 



0,1 

0,2 

0,3 

0,4 

0,5 


0,1 

0,230274 

0,160403 

0,109667 

0,071529 

0,043153 


0,2 

0,165211 

0,103981 

0,071011 

0,042870 

0,023259 

1,0 

0,3 

0,119886 

0,075975 

0,045964 

0,025566 

0,012414 


0,4 

0,086363 

0,051617 

0,029083 

0,014798 

0,006389 


0,5 

0,060786 

0,033833 

0,017533 

0,008069 

0,003067 


0,1 

0,294393 

0,202735 

0,137144 

0,089016 

0,053983 


0,2 

0,226947 

0,149552 

0,095926 

0,058390 

0,032695 

1,5 

0,3 

0,174787 

0,109232 

0,065923 

0,037363 

0,019156 


0,4 

0,132053 

0,077255 

0,043329 

0,022596 

0,010496 


0,5 

0,096059 

0,051488 

0,026291 

0,012369 

0,005106 


0,1 

0,312157 

0,206336 

0,134101 

0,084080 

0,049815 


0,2 

0,245598 

0,154311 

0,094757 

0,055747 

0,030645 

2,0 

0,3 

0,191602 

0,113072 

0,064891 

0,035454 

0,017914 


0,4 

0,145469 

0,079205 

0,041741 

0,020831 

0,009516 


0,5 

0,105171 

0,051312 

0,024126 

0,010695 

0,004302 


0,1 

0,306503 

0,191072 

0,117287 

0,069869 

0,039687 


0,2 

0,241768 

0,142001 

0,081812 

0,045535 

0,023987 

2,5 

0,3 

0,188126 

0,102551 

0,054663 

0,028048 

0,013516 


0,4 

0,141580 

0,070018 

0,033780 

0,015622 

0,006767 


0,5 

0,100566 

0,043493 

0,018294 

0,007339 

0,002723 


0,1 

0,290584 

0,168953 

0,096914 

0,054251 

0,029211 


0,2 

0,228110 

0,123742 

0,066078 

0,035928 

0,017053 

3,0 

0,3 

0,175929 

0,087443 

0,042699 

0,020240 

0,009152 


0,4 

0,130526 

0,057832 

0,025128 

0,010600 

0,004283 


0,5 

0,090664 

0,034259 

0,012691 

0,004476 

0,001014 
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